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Abstract 

ISO has enabled far-infrared spectroscopy of a vari- 
ety of galaxies. Using the [CII] (158/xm) and [OI](63/xm) 
lines, we can characterize the physical conditions in the 
star-forming ISM. These observations also form the basis 
of our expectations for what the Herschel Space Obser- 
vatory will observe for high redshift galaxies. While [CII] 
is suppressed in ULIRGs and normal galaxies with high 
dust temperatures, it is stronger than expected in metal 
poor galaxies by factors of a few. Young galaxies at high 
redshifts might be expected to be both metal poor and 
actively star-forming, leading to contrary expectations for 
the [CII] line strength. The best prediction for [CII] de- 
tection is derived by using the observed proportionality 
between [CII] and mid-IR emission from PAHs. Using the 
observed [CII]/7 /xm ratio and number counts from ISO 
deep surveys we predict that HSO will be able to detect 
100 sources/square-degree in the [CII] line. 
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1. Introduction 

Herschel Space Observatory (HSO) will have the wave- 
length coverage and sensitivity to detect high redshift galax- 
ies at the peak of their dust continuum emission. The three 
instruments onboard HSO will also carry out spectroscopy 
in the far-infrared and sub-millimeter. From the FIR lines 
of C + and O, we can derive physical conditions in the 
star-forming ISM at these redshifts. From optical obser- 
vations it seems that star-formation peaks at z=l-2 (e.g. 
Madau et al. 1996). It would be interesting to investigate 
the properties of the ISM in star-forming galaxies at these 
redshifts, and thus learn the causes and consequences of 
higher star-formation. 

In this paper I use ISO observations of the fine struc- 
ture lines in nearby galaxies to (a) summarize what we 
have learned about the ISM of galaxies from ISO and (b) 
do feasibility calculations about what HSO will see at high 
redshifts. 



2. The [CII] (158/im) line 

[CII] (158/xm) line is the dominant coolant of the neutral 
ISM in all but the hottest galaxies. For most of the ob- 
served galaxies 0.1-1% of the FIR continuum emerges in 
this single line. Not coincidentally, it is the best studied 
line in this wavelength regime. It has been observed for 
normal galaxies (Stacey et al. 1991, Malhotra et al. 1997, 
1999, 2001, Leech et al. 1999, Pierini et al. 1999), irregu- 
lar/dwarfs (Smith et al. 1997, Bergvall et al. 2000, Hunter 
et al. 2001, Madden et al. 2001), Ellipticals (Malhotra et 
al. 2000, Unger et al. 2000) and luminous and ultralumi- 
nous infra-red galaxies (Luhman et al. 1998, Fischer et al. 
2001). So now we have a low redshift "basis set" of various 
types of galaxies to draw upon. 

The deficiency in [CII] /FIR in ULIRGs and normal 
galaxies (Malhotra et al. 1997, Luhman et al. 1998) came 
as a surprise to many. The decrease in [CII] /FIR correlates 
best with the IRAS colors of galaxies F zy (60/im)/F zy (100/xm) 
(Figure 1; Malhotra et al. 2001). It also correlates, but less 
strongly, with other quantities like FIR/Blue colors and 
Infrared luminosity. Since F^(60/im)/F t/ (100//m), FIR/Blue 
and FIR luminosity correlate with each other, we cannot 
say which correlation with [CII] /FIR is primary and which 
ones are secondary. Luminous and ultraluminous galaxies 
from the sample of Luhman et al. 1998 follow the same 
trends. Two high redshift quasars BRI 1202-0752 (z=4.69) 
and BRI 1335-0415 (z=4.41) have measured upper limits 
on the [CII] flux and also follow the same trendd (Bcnford 
1999). 



3. Irregular galaxies 

While FIR colors F„(60/xm)/Fi,(100/xm) show the strongest 
correlation with L[ch]/Lfir, there seems to be a second 
parameter, which is apparent from Figure 1. Irregular 
galaxies have a higher Lrcni/LprR and Ellipticals have 
lower Lrcni/LpiR ratio. It is not yet clear whether this 
has to do with the lower metallicity affecting the chem- 
istry of the ISM directly or because low metallicity stellar 
populations produce a harder radiation field. Hunter et al. 
2001 and Madden et al. 2001 discuss in detail how both 
L[ch]/Lfir and [CII] /MIR emission is higher for irregular 
galaxies. 

Even more dramatic is the fact that [OIII] (88/xm) line 
from HII regions is very bright in 2 of our irregular galax- 
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Figure 1. The ratio of [CII] to far-infrared continuum, 
L[ C h]/Lfir, * s plotted against the ratio of flux in the IRAS 
60/xm and lOOimr bands, F„(60iim)/F„(100/xm). Filled 
circles are normal galaxies from the ISO-Key Project sam- 
ple. Irregular galaxies are denoted with the star sign and 
ellipticals with bulls eye symbols. Triangles are luminous 
and ultraluminous galaxies from the sample of Luhman 
et al. (1998). There is a trend for galaxies with higher 
F t/ (60/im)/F^(100/im) (indicating warmer dust) to have 
lower L[ch]/Lfir, for normal as well as ULIRGs. Two 
normal galaxies in a sample of 60 have no detected [CII], 
and they are identified with labels and shown as upper 
limit symbols; other upper limits come from Luhman et 
al. (1998). Rank correlation tests show that Lpjjj/LpiR 
and Fi / (60/im)/Fi / (100/im) are correlated at the 4-4°~ level. 
L[cii]/Lfir upper limits for two high redshift sources also 
follow this trend (Benford 1999). 



ies. In IC 4662 and NGC 1569 [OIII]/FIR is fully 1% (Mal- 
hotra et al. 2001). This is promising for detection of [OIII] 
at high redshifts, but it also means that we cannot reliably 
assume that that the brightest (and sometimes the only) 
line detected is [CII]. 

NGC 1569 and IC 4662 are dwarf galaxies and hence 
not very luminous. We would not be able to see thier coun- 
terparts at z=l or even z=0.5. The [OIII] (88/j,m) line from 
the metal poor irregular Haro 11 (Bergvall et al. 2000) 
would be observable at z=l, even as the [CII] line falls 
short. 



4. FIR SPECTROSCOPY AT HIGH REDSHIFTS 

In figure 2 we compare the luminosity distribution of three 
prominent FIR lines in the ISO key-project sample. These 
lines are [CII] (158 /jm), [OI] (63 /im) and [OIII] (88 /mi). 
This is not a luminosity function, since our sample is not 
volume limited, but serves to compare the luminosity dis- 
tribution among the three brightest lines. We see that in 
spite of the [CII] deficiency which cuts off the high lu- 
minosity tail of the distribution, [CII] is clearly the more 
luminous of these three lines, and therefore offers the best 
prospects for being widely detected. 
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Figure 2. The luminosity distribution of the three promi- 
nent FIR lines [CII] (158 fim), [OI] (63 /mi) and [OIII] 
(88 /im,). The sample here is the ISO Key-Project on nor- 
mal galaxies. Since this sample was designed to span the 
possible galaxy properties rather than be a representative 
or a volume limited sample, this figure does not show the 
luminosity function in any of the lines. However it is in- 
structive to see the relative luminosity distribution in the 
lines and to note that [CII] is still the most luminous and 
likely to be detected at redshifts 0.5-1. The x-axis shows 
the log of the fluxes (in units of 1 x lQ~ 20 W/m 2 ) of these 
galaxies if they were at a distance of 1 Gpc 

The [CII] deficiency in the most active and luminous 
galaxies flattens the luminosity function of the line rela- 
tive to continuum. Since the most luminous objects (e.g. 
ULIRGs) are deficient in [CII] the prospects for detect- 
ing them at redshifts « 4 and higher are dim. Observa- 
tions of known high redshift objects with JCMT and CSO 
have only yielded upper limits (van der Werf et al. 1998, 
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Bcnford 1999). These upper limits are consistent with the 
trend seen between L[ch]/Lfir and F„(60/im)/F„(100^tm) 
(Figure 1). It seems that the correlation between L[ch]/Lfir 
and F 1/ (60^m)/F 1/ (100/im) is stronger than between L[cn]/L 
and luminosity (Malhotra et al. 2001). In principle, a cool 
but luminous source would be detectable at high rcdshifts. 
In practice, however, luminous sources also tend to have 
warm FIR colors. The [CII] deficiency in high redshift 
galaxies may be somewhat mitigated if starburst galaxies 
are metal poor, since metal poor galaxies have [CII] /FIR 
higher by a factor of a few. But this effect is not expected 
to be large. 

4.1. Source density of high-z [CII] targets 

While the [CII] /FIR behaviour is complicated and corre- 
lates best with F„(60/zm)/F„(100/mi), we can exploit the 
[CII] /MIR constancy to predict the number of galaxies 
that are observable with Herschel Space Observatory. 

Hclou ct al. 2001 point out that the ratio of [CII] flux 
and the flux in the mid-infrared band at 7 yum remains 
constant with FIR colors and shows a smaller scatter than 
does Ljchj/Lfir. The physical interpretation of the cor- 
relation between [CII] and 7/j.m flux is simple. The mid- 
infrared flux in the 6.75 band of ISO-CAM is dominated by 
emission from Polycyclic Aromatic Hydrocarbons (PAHs). 
We interpret the stable [CII]/F 7Mm ratio as evidence that 
gas heating is dominated by the PAHs or small grains 
which are also AFE carriers. 

Another interpretation of the constancy of [CII]/F 7Atm 
is that PAHs/ very small grains and C + are co-extensive in 
PDR regions. The decrease in both [CII] /FIR (Malhotra 
et al. 1997, 2001) and F 7Aim /FIR with F60/F100 (Helou et 
al. 2001, Lu et al. 2001) are then due to a smaller fraction 
of FIR arising from the PDR phase. 

Regardless of the right explanation/interpretation of 
the proportionality between [CII] and PAH emission, we 
can use it to predict what the Herschel Space Observatory 
will be able to observe. 

From the ISO key-project sample we derive 

[CII] = (Fr^/lmJy) x W- 17 - 2 W/m 2 

A sensitivity of 6 x IQ-^W/m 2 (HIFI) in the [CII] line 
then translates to rest-frame 7 micron flux of 1 mJy. Deep 
observations with ISOCAM (Elbaz et al. 1999) show that 
the surface density of sources with 15 ^m fluxes > lmJy 
is roughly 100 per square-degree. 

There are two caveats we need to keep in mind when 
using the deep counts in the mid-IR to predict the num- 
ber of galaxies which will have detectable [CII] flux. The 
first is that the flux in the mid-infrared is dominated by 
spectral features due to aromatic molecules, which means 
that as we get to higher rcdshifts K-corrections for fluxes 
measured in any broad-band filter are substantial and not 
always of the same sign (Xu et al. 1998). The other con- 
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Figure 3. The ratio of [CII] flux and 7 /im flux is un- 
varying with FIR colors and shows a smaller scatter than 
L [C ii]/Lfir (a([CII\/F riaa ) = 0.22. The [CII] flux is m 
units of W~ 14 W/m 2 and the 7 /im flux is in mJy. 



cern has to do with metal poor galaxies. They show a 
deficiency in PAH features in the mid-IR (Thuan et al. 
1999), but also show higher [CII] /MIR and [CII] /FIR by 
factors of 2-3. In these galaxies presumably most of the 
heating is due to the very small grain component of dust, 
which emit in the mid-IR but do not produce aromatic 
features. 

5. Physical conditions in the PDRs 

We derive the average physical conditions in the neutral 
gas in galaxies, by comparing the observed line to contin- 
uum ratio (L[cn] +L[oi])/Lfir and F !/ (60^m)/Fj / (100/im) 
with PDR models (e.g. Kaufman et al. 1999). The main 
results of this study are (Malhotra et al. 2001): 

(1) The derived temperatures at the PDR surfaces 
range from 270-900 K, and the pressures range from 6 x 
10 4 — 1.5 x 10 7 Kcm~ 3 . The lower value of the pressure 
range is roughly twice the local solar neighborhood value 
and the upper end is comparable to pressures in HII re- 
gions in starburst galaxies (Heckman, Armus & Miley 
1990) which also corresponds to the pressure and sur- 
face brightness at which starbursts saturate (Meurer et 
al. 1997). 

(2) The average FUV flux Go and gas density n scale 
as Go oc n 1A . This correlation is most naturally explained 
as arising from Stromgrcn sphere scalings if much of the 
line and continuum luminosity arises near star-forming 
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regions. From simple Stromgren sphere calculations (cf. 
Spitzer 1978) we can derive that the FUV flux at the neu- 
tral surface just outside the Stromgren sphere should scale 
as Go oc n 4 / 3 , which is consistent (within errors) to the 
scaling seen in Figure 4. The Go, n and P which we derive 
for given galaxy represent a luminosity weighted average 
value and mostly represents dense GMCs which lie close 
to the OB stars. 

(3) The range of Go, n and P which we derive for the 
different galaxies can reflect several interesting differences 
in their star formation processes and histories. If global 
star formation is episodic then high Go and n imply that 
the galaxy is observed shortly after a burst because the 
OB stars have not moved far from their natal clouds. Al- 
ternatively, the differences in Go, n and P may reflect 
differences from galaxy to galaxy in the properties of the 
GMCs which form the OB stars. Larger GMCs may keep 
their OB stars embedded for a longer fraction of their life- 
time, resulting in higher average Go and n. Or the GMCs 
could be the same size but denser, so that the higher den- 
sity ambient gas would lead to smaller Stromgren spheres 
and higher Go at the edges of the spheres. 



6. Conclusions 

It seems unlikely that spectroscopy with HSO will be used 
to find redshifts of infrared bright galaxies with unknown 
or obscured optical counterparts simply because scanning 
the possible redshift range would be time consuming. But 
once the redshifts are known, one can use the observations 
of [CII] and [01] (63 /jm) lines to derive the physical con- 
ditions in PDRs. Or if photometric redshifts using mid to 
far-infrared colors can be refined enough, one could search 
for fine structure lines to get the precise redshifts. 

FIR spectroscopy of modest redshift ( z — 1) galaxies 
would be feasible with the HSO. This would be valuable 
to understand the higher global rates of star-formation in 
this epoch at roughly half the age of the universe. With 
ISO we have been able to characterize the ISM in star- 
forming local galaxies, so we do have a comparison set. 
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Figure 4- This figure shows the derived far- UV flux Go and 
gas density n solution for the ISO key-project galaxies 
based on comparison of FIR line data and PDR models 
of Kaufman et al. 2001. A least square fit is made to the 
Go vs. n relation assuming equal error in both axes. The 
best fit slope is 1.4, i.e. Go oc n 1A , which is consistent 
with PDRs surrounding ionization bounded expanding HII 
regions. 
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